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Abstract. We consider the two-body problem in post-Newtonian approximations of gen- 
, eral relativity. We report the recent results concerning the equations of motion, 

O" 1 and the associated Lagrangian formulation, of compact binary systems, at the 

»' | | third post- Newtonian order (~ 1/c 6 beyond the Newtonian acceleration). These 

(3J[), equations are necessary when constructing the theoretical templates for searching 

and analyzing the gravitational-wave signals from inspiralling compact binaries in 
VIRGO-type experiments. © 2001 Academie des sciences/Editions scientifiques 
et medicales Elsevier SAS 



general relativity /equations of motion 

Sur le probleme des deux corps en relativite generate 

Resume. Nous considerons le probleme des deux corps dans V approximation post-newtonien- 
ne de la relativite generate. Nous presentons les resultats recents concernant les 
equations du mouvement, et la formulation lagrangienne associee, de systemes bi- 
naires d'objets compacts, au troisieme ordre post-newtonien 1/c 6 apres I' acce- 
leration newtonienne) . Ces equations sont utilisees dans la construction des fil- 
tres theoriques pour la detection et I 'analyse des signaux d'ondes gravitationnelles 
provenant des binaires compactes spiralantes dans les experiences du type VIRGO. 
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By two-body problem we mean the problem of the dynamics of two strutureless, non-spinning 
point-particles, characterized by solely two mass parameters mi and m2, moving under their 
mutual, purely gravitational, interaction. Surely this problem, because of its conceptual simplicity, 
is among the most interesting ones to be solved within any theory of gravity. Actually, there are 
two aspects of the problem: the first sub-problem consists of obtaining the equation of the binary 
motion, the second is to find the (hopefully exact) solution of that equation. We refer to the 
equation of motion as the explicit expression of the acceleration of each of the particles in terms 
of their positions and velocities. It is well known that in Newtonian gravity, the first of these 
sub-problems is trivial, as one can easily write down the equation of motion for a system of N 
particles, while the second one is difficult, except in the two-body case N = 2, which represents in 
fact the only situation amenable to an exact treatment of the solution. 

In general relativity, even writing down the equations of motion in the simplest case N = 2 is 
difficult. Unlike in Newton's theory, it is impossible to express the acceleration by means of the 
positions and velocities, in a way which would be valid within the "exact" theory. Therefore we are 
obliged to resort to approximation methods. Let us feel reassured that plaguing the exact theory 
of general relativity with approximation methods is not a shame. It is fair to say that many of 
the great successes of this theory, when confronted to experiments and observations, have been 
obtained thanks to approximation methods. Furthermore, the beautiful internal wheels of general 
relativity also show up when using approximation methods, which often deserve some theoretical 
interest in their own, as they require interesting mathematical techniques. 

Here we whall investigate the equation of the binary motion in the post-Newtonian approxima- 
tion, i.e. as a formal expansion when the velocity of light c tends to infinity, up to the so-called 
third post-Newtonian (3PN) order, i.e. 1/c 6 beyond the usual Newtonian acceleration. This prob- 
lem is not pure academicism, because the equation of motion at this somewhat frightening 3PN 
order is needed (among other things) during the computation of the orbital phase of the inspi- 
ralling compact binaries [JjJ |). The phase constitutes the crucial component of the theoretical 
templates to be used for searching and analyzing the gravitational waves from the binary inspiral 
in VIRGO-type detectors. 

We write the 3PN acceleration of the particle 1, say, in the form 



a l — + ^2^1PN + ~[ ^2PN + ^5^2.5PN + ^3 ^3PN + ® ^^7 ) 1 1 1 

where the first term is given by the famous Newtonian law, 
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For simplicity, we do not write the particle's label 1 on the A^pn's. We denote by r i2 = |yi(i) — 
y2 (£) | the distance between the two particles, with yi = (y\) and y 2 = (y|) their "absolute" 
positions, by n\ 2 = (y\ — V % 2)l T n the corresponding unit direction, and by a\ = dv\/dt and 
v\ = dy\/dt the acceleration and velocity of 1 (and idem for 2). Sometimes we pose v\ 2 — v\ — v 2 
for the relative velocity; and similarly for the relative position and acceleration. The equation for 
the other body is obtained by exchanging the labels 1 <-> 2 (remembering that n\ 2 changes sign in 
this operation). 

As a consequence of the equivalence principle, which is incorporated "by hand" in Newton's 
theory and constitutes the fundamental basis of general relativity, the acceleration of 1 should 
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not depend on mi (nor on its internal structure), in the "test-mass" limit where the mass mi is 
much smaller than mi- This is of course satisfied by the Newtonian acceleration (||), which is 
independent of mi, but this leaves the possibility that the acceleration of the particle 1, in higher 
approximations, does depend on mi, via the so-called self- forces, which vanish in the test-mass 
limit. Indeed, this is what happens in the post-Newtonian corrections computed in Eqs. (|^)-@, 
which show explicitly many terms proportional to (powers of) m\. 

Though the successive post-Newtonian approximations are really a consequence of general rel- 
ativity, they should be interpreted using the common-sense language of Newton. That is, having 
chosen a convenient general-relativistic (Cartesian) coordinate system, like the harmonic coordi- 
nate system adopted below, we should express the results in terms of the coordinate positions, 
velocities and accelerations of the bodies. Then, the trajectories of the particles can be viewed as 
taking place in the absolute Euclidean space of Newton, and their (coordinate) velocities as being 
defined with respect to absolute time. Not only this interpretation is the most satisfactory one 
from a conceptual point of view, but it represents also the most convenient path for comparing the 
theoretical predictions and the observations. For instance, the solar system dynamics at the first 
post-Newtonain level is defined, following a recent resolution of the International Astronomical 
Union, in a harmonic coordinate system, the Geocentric Reference System (GRS), with respect to 
which one considers the "absolute" motion of the planets and satellites. 

With the same spirit we present below the 3PN equations of motion in terms of the harmonic- 
coordinate positions and velocities of the particles, i.e., in particular, in a form which is not 
manifestly covariant. But because the equations come from general relativity, they are endowed 
with the following properties, which make them truly "relativistic" . 

(I) The one-body problem in general relativity corresponds to the Schwarzschild solution, so the 
equations possess the correct "perturbative" limit, that given by the geodesies of the Schwarzschild 
metric (developed to the 3PN order), when the mass of one of the bodies tends to zero. 

(II) Because general relativity admits the Poincare group as a global symmetry (in the case of 
asymptotically flat space-times), the harmonic-coordinate equations of motion stay invariant when 
we perform a global Lorentz transformation, expanded at the 3PN order. 

(III) Since the particles emit gravitational radiation there are some terms in the equations which 
are associated with radiation reaction. These terms appear at the order 2.5PN or 1/c 5 (where 
5 = 2s + 1, s = 2 being the helicity of the graviton); see Eq. ([!]). They correspond to an "odd"- 
order post-Newtonian correction, which does not stay invariant in a time reversal. By contrast, as 
we shall see, the "even" -orders 1PN, 2PN and 3PN correspond to a dynamics which is conservative. 

(IV) General relativity is a non linear theory (even in vacuum) , and some part of the gravitational 
radiation which was emitted by the particles in the past scatters off the static gravitational field 
generated by the rest-masses of the particles, or interacts gravitationally with itself. The "tail" 
radiation, produced by non-linear scattering, converges back onto the system at our present epoch, 
and modifies its current dynamics. The effect appears at the 4PN order [negligible in Eq. ([!])], 
and corresponds to some "even" -order modification of the radiation-reaction damping j3|. 

The dominant relativistic correction term (1PN, or 1 /c 2 ) in the case of two non-spinning compact 
bodies was derived first by Lorentz and Droste Subsequently, Einstein, Infeld and Hoffmann 
j^] obtained the 1PN corrections by means of their famous "surface-integral" method, in which 
the equations of motion are deduced from the vacuum field equations. This method is applicable 
to any compact objects (they be neutron stars, black holes, or, maybe, naked singularities). The 
lPN-accurate equations were also obtained, for the motion of the centers of mass of extended 
bodies, by Fock and followers using a technique that can be qualified as more "physical" than 
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the surface-integral method, as it takes explicitly into account the internal structure of the bodies. 
The result is 
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A(n 12 vi) - Z(n 12 v 2 ) 
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We denote the usual Euclidean scalar product by parenthesis, for instance (ni 2 vi) — ni2.vi and 
(viv 2 ) = V1.V2. Witness the first term in Eq. (|J), proportional to mi, which represents a self- force 
at the 1PN order. The lPN-accurate equations of motion (known also for N compact bodies) are 
sometimes called after the names of Einstein, Infeld and Hoffmann || . 

Things become more and more difficult when going to higher and higher post-Newtonian ap- 
proximations. A technical problem is the rapid proliferation of terms. Typically any allowed term 
(compatible dimension, correct mass dependence) does appear with a definite non-zero coefficient 
in front. As a result, the expressions look quite complicated; that's maybe the price we have to pay 
for expressing in a Newtonian fashion what really is a relativistic equation. But, the point for our 
purpose is that the Newtonian- like equations are fully explicit: the successive relativistic correc- 
tions depend only on the instantaneous positions and velocities — all accelerations, appearing at 
intermediate stages of the calculation, are consistently "order-reduced" by means of the equations 
of motion themselves. 

The 2PN approximation was tackled by Otha et al [Q, how considered the post-Newtonian 
iteration of the Hamiltonian of N point-particles. We refer here to the Hamiltonian as the Fokker- 
type Hamiltonian, which is obtained from the matter-plus-field Arnowitt-Deser-Misner (ADM) 
Hamiltonian by eliminating the field degrees of freedom. The result for the 2PN equations of 
binary motion in harmonic coordinates was obtained by Damour and Deruelle |8|, |l0| , building 
on a non- linear iteration of the metric of two particles initiated in Ref. fill]. It reads 
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The corresponding result for the ADM-Hamiltonian of two particles at the 2PN order was given 
by Damour and Schafer [12]. 
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As mentioned above, the equation of motion is conservative up to the 2PN level. Only at the 
2.5PN order appears the first non-conservative effect, associated with the gravitational radiation 
reaction. The (harmonic-coordinate) equation of motion up to that level has been derived by 
Damour and Deruelle ^ ^, [l(| , who used it for the study of the radiation damping of the orbital 
period of the binary pulsar PSR 1913+16 [^0|. In addition, the 2.5PN-accurate ADM-Hamiltonian 
was worked out by Schafer |l5)| , and the harmonic-coordinate equations as well as the complete 
gravitational field of two particles were computed by Blanchet, Faye and Ponsot |l4j by means of a 
new method coming from the work on gravitational- wave generation jj| . The latter calculations use 
a formal description of the objects by delta functions. Needless to say, there is mutual agreement 
between all the results obtained so far. The 2.5PN radiation damping term reads 
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It is important to realize that the 2.5PN equations of motion have been proved to hold in 

the case of binary systems of strongly self-gravitating bodies ]l(J . This is via an effacing principle of 
the internal structure of the bodies. As a result, the equations depend only on the "Schwarzschild" 
masses, m\ and m 2 , of the compact objects. Compacity parameters, namely the ratios and 

j^j- between the masses and the radii b\ and b 2 of the objects, do not enter the equations of 
motion. This fact has been explicitly verified at the 2.5PN order by Grishchuk and Kopejkin |15| ], 
who have made a "physical" computation, a la Fock, taking into account the internal structure 
of two self- gravitating extended bodies. Furthermore, the same 2.5PN equations of motion have 
also been established by Itoh, Futamase and Asada |ll| , who use a variant of the surface- integral 
approach g], valid for compact bodies, independently of the strength of the internal gravity. 

The present state of the art is the 3PN approximation. The equations to this order have been 
worked out independently by two groups, by means of different methods, and with equivalent 
results. On one hand, Jaranowski and Schafer fll7||, and Damour, Jaranowski and Schafer JlT 
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employ the ADM-Hamiltonian formalism of 



general relativity; on the other hand, Blanchet and Faye |20|, |2lJ], and Andrade, Blanchet and 
Faye J2^] , founding their approach on the post-Newtonian iteration initiated in Ref . Hj] , compute 
directly the harmonic-coordinate equations of motion. The end results are physically equivalent 
in the sense that there exists a unique "contact" transformation of the dynamical variables, that 
changes the harmonic-coordinate Lagrangian obtained in Ref. |^2| into a new Lagrangian, whose 
associated Hamiltonian coincides exactly with the one given in Ref. |iiq] . The 3PN contribution 
to the acceleration of the particle 1 reads, in a Cartesian harmonic coordinate system, 
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Notice that the 3PN term involves some 7r 2 . Typically, the ir 2 terms arise from non-linear interac- 
tions involving some integrals such as — f -i 2 ? = — . 

More importantly, we see that the 3PN acceleration depends on three arbitrary constants: two 
length scales r[ and r' 2 entering some logarithms, and a dimensionless constant A. It was proved 
in Ref. ]^0) that r[ and r 2 are merely linked with the choice of harmonic coordinates. Indeed, as 
we are using point-particles, the usual condition of harmonic coordinates, viz d v W — 0, does not 
completely fix the gauge, as we can always add a gauge vector £ M — 5x^ , which satisfies A£ M = 
and is singular at the location of the two particles. The constants r[ and r' 2 can thus be removed 
by a coordinate transformation; as a result they will never appear in any physical result, such as 
the invariant center-of-mass energy of the binary. 

By contrast with the harmless gauge-constants r[ and r' 2 , the constant A represents a true phys- 
ical ambiguity, which reflects probably an incompleteness of the Hadamard regularization used to 
cope with the infinite self-field of point-particles. This regularization is based on Hadamard's con- 
cept of the "partie finie" of singular functions and divergent integrals. Actually, it has been found 
necessary for solving this problem to develop an extended version of the Hadamard regularization, 
and a theory of generalized functions associated with it [Ml . Very likely, the presence of A is related 
to the fact that, starting from the 3PN order, many non-linear integrals composing the equations 
of motion, when taken separately, depend on the internal structure of each body, even in the limit 
where the radius tends to zero. The results given by the ADM-Hamiltonian approach Jlj], IS] 
depend also on one arbitrary physical parameter, called w sta tic- More precisely, the authors [17] 
introduced originally two unknown constants, w s tatic and ^kinetic, but akinetic was fixed later to a 
unique value by imposing, in an ad hoc manner, the global Poincare invariance of the Hamiltonian 
flS[| . On the other hand, the authors 20, 2l| have only one constant A because their regularization 
is defined in a Lorentz-invariant way. The equivalence between the approaches |2^, |2l], |2^] and 
fl7| , |l8f holds if and only if A = — fj^static — 'gif^- Recently, the value ^static = has been obtained 
by means of a different regularization (dimensional) within the ADM-Hamiltonian approach IljJ]. 
This result would mean that A = — (but we keep A unspecified in the present discussion). 

Going to still higher post-Newtonian orders, we mention that the subdominant radiation-reaction 
effect at the 3.5PN order has been calculated in Ref. |2j| for two bodies in an arbitrary gauge, and 
in Ref. ]^4| for general "fluids" in a Burke-Thorne-extended gauge (there is agreement between 
the two methods (2j|). Also, we know the contribution of the tails of waves in the equations of 
motion, which appears at the 4PN, or 1/c 8 , order ||. 

The complicated expressions (0)-(||) simplify drastically in the case of an orbit which is circular, 
apart from the gradual inspiral driven by radiation reaction. This case corresponds to the physical 
situation of most inspiralling compact binaries, since the radiation reaction forces tend to circularize 
rapidly the orbit. In this case the relative acceleration reads 
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Relative quantities are denoted e.g. by y\ 2 — y\ — y 2 ; mass parameters are the total mass m — 
mi + m2 and the mass ratio v = mim 2 /m 2 such that < v < |, with v = j when the two 
masses are equal, and v — > in the test-mass limit for one of the particles. The second term in 
Eq. (0) is the radiation reaction force, whose expression follows immediately from Eq. (||), while 
to represents the orbital frequency of the circular motion at the 3PN order, and is related to the 
harmonic-coordinate separation r\ 2 by the "Kepler" law 
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The post- Newtonian corrections are parametrized by 7 = r ^™y ■ The constant r' Q is given by In r' = 
^\nr[ + 22- In 7*2, where and r' 2 are the two gauge-constants in Eq. (|^). The circular-orbit 
equation (0)-(|^) is used, notably, to compute the third time-derivative of the binary's quadrupolc 
moment, which constitutes the main contribution to the gravitational-wave flux at infinity, at the 
3PN order §. 

Finally, coming back to the general case of non-circular orbits, we show that the 3PN equations 
of motion, when neglecting the radiation-reaction term at the 2.5PN order, admit a Euler-Lagrange 
formulation. As it turns out, the Lagrangian in harmonic coordinates is a generalized one, in the 
sense that it depends not only on the positions and velocities of the particles, like an ordinary La- 
grangian, but also, starting at the 2PN order, on their accelerations. This fact has been established 
in Refs. H [n], [l2|]. Thus, the equations of motion take the form 

m x 4. 5PN = o(i?) > (9) 

where A l 2 5PN is given by Eq. (^|), and where the 3PN harmonic-coordinate generalized Lagrangian 
is of the type 

L = Ln H 2-^lPN + ~ ^2PN + — L 3 pn + O ( — J . (10) 

The Newtonian piece reads 

mivf m 2 v\ Gm 1 m 2 . . 

= — 7, 1 o 1 • l n J 

2 2 ria 

The 1PN piece, which is still "ordinary" (not depending on accelerations), has been obtained by 
Lorentz and Droste 101, and Fichtcnholz [E5|: 
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To the expression given above, one must add the terms corresponding to the label exchange 1 «-> 2, 
including those that are already symmetric under exchange. Now, starting at the 2PN order, we 
get a dependence over the accelerations || [h], ^2|. By adding to the Lagrangian some so-called 
"multi-zero" terms, which do not contribute to the equations of motion, one can always arrange 
that the dependence over the accelerations be linear. Of course, it is not allowed to replace the 
accelerations by the equations of motion in a Lagrangian (however, this can and should be done 
in the final expressions of the conserved integrals derived from that Lagrangian). We get 
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Finally, the 3PN-accurate piece of the harmonic-coordinate Lagrangian depends also on accelera- 
tions; it is notable that accelerations are sufficient, there is no need to include, at the 3PN order, 
derivatives of accelerations. We find p2] | 

5 o G 2 m\m 2 (Yi 4 83. .0. . 35. , 2 . , 2 

^3PN = Tr>o m l V l ^ 2 T^("l2^l) + 7^(ni2«l) \ni2V2) - -^-(ni2«l) (ni2^2j 

izo r 12 \to to 

245, , 2 2 179. ,, . 2 235. , 2 2 373 4 529. , 2 . 
- -^{n 12 vi) v x + —{ni2Vi){rii2V 2 )v 1 - —{n 12 v 2 ) v x + — v x + — (m 2 Ui) («iUa) 

97 719 463 7 

y(ni2fi)(ni 2 v 2 )0iV2) ~ -^l^a) + -^-{viv 2 ) 2 - — (n 12Vl ) 2 

1, . o 1, ,0 o 463 , , 19, , o 45 , 



1 , \ t X 2 1, , 22 « 2 2 19, , 2 45 4 \ 

- 2<. ni2i;i K ni2W2 ) u 2 + ^(ni2«2j v 2 + ^" v i v 2 " y(«i«2)w 2 + ^2 I 

/3 5 1 

+ Gm 1 m 2 f -(a 1 u 2 )(ni2fi)(ni2W2) 2 + j^(aiv 2 )(n 12 v 2 ) 3 + -(n 12 a 1 )(n 12 Ui)(ni2W2) 3 

+ -7:(ni2ai)(«i 2 W2) 4 + ^-(aiv 1 )(n 12 v 2 )v' 2 - (ai« 2 )(rii2f2)«i - 2(a 1 v 1 )(n 12 v 2 )(v 1 v 2 ) 
16 4 

13 5 

+ -j(aiv 2 )(ni 2 v 2 )(viv 2 ) + -(n 12 ai){ni 2 v 2 ) 2 (viv 2 ) - -{ni 2 ai){ni 2 vi) 2 vl 

4 o o 

15 15 1 

+ — {aivi){ni 2 v 2 )vl - —{aiv 2 )(ni 2 v 2 )vl - -(ni 2 ai)(ni 2 Vi)(ni 2 v 2 )v 2 
o o Z 

- ^(«i2ai)(ni 2 W2) 2 «|^ + G ™ lJn2 ( - ^(a 2 Ui)(7ii2Ui) - ^{ni 2 a 2 ){ni 2 vi) 2 

~ -7rr( a i v 2){ni2V 2 ) - ^r("-i 2 ai)(ni2W 2 ) 2 + ^r(ni2ai)«i - ^-{ni 2 a 2 )vf 
24 6 16 48 

185. 20. . 2 \ Gm 1 m 2 ( 5 , 

— (ni2aij(uiw 2 ) + — (ni 2 ai)w 2 M — (n 12 ui) (ni 2 v 2 ) 

o 6 J ri2 \ 0/ 

\{n 12 vi){ni 2 v 2 )' i vl + ^(n 12 v 2 ) 4 v 2 - ^-{ni 2 vi){ni 2 v 2 )v\ + ^-(n 12 v 2 ) 2 v* + ^vf 

5 o lb 4 lb 

L5. .0. >o. w \ o / \ 



8 V M ; 1 8 V ; 1 16 v M 
15 

— («i2«i) 2 (^i2f2) 2 (vif2) + {ni 2 vi)(ni 2 v 2 )v 2 (viv 2 ) - 

13 5 1 

—vf(viv 2 ) + —{n 12 vi)(n 12 v 2 )(viv 2 ) 2 + — (wiw 2 ) 3 
16 16 16 

P(ni2fi)(ni2W 2 )w 2 « 2 + ^1^2 - j^l^l^)"^ 

3G 4 mfm 2 G A ra\m\ ( 5809 11 22 /r. 



ri2 V 

'(ni 2 i;i)(ni2«2)«i + ^("-i2«2) 2 
3 

-z{nx 2 v 2 ) 2 v\{vxv 2 ) 

o 

^{n 12 v Y ) 2 v\vl 
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G 3 mfml f383^ 2 889 305 2 439 
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— (nia«i) - -^(ni2V 1 )(n 1 2V 2 ) - -^■v 1 + — («iu 2 ) 



123 , >9 n 123 , . , \ 2 n 9 ^""^ / \ 2 

— (ni2« 1 ) 7T + — (ni2Wi)(ni2V 2 )7r + — ^tt - — (v^tt 
G 3 m\m 2 ( 8243^ 2 15541 3 2 



12 



21Q (niaui) + -^2Q-(ni2Vi)(ni2U2) + ^O 1 ^) 



15611 2 17501 , N , 5 2 , „„, w N1 fr 12 



+ (v-\V2) -\ — w 9 + 22(ni 2 Wi )(ni2Vi2) In , 

1260 1 1260 v ; 4 2 v 12 in 12 12; , 



See Ref. |22 for the explicit expressions of the ten conserved quantities, at the 3PN order, corre- 



sponding to the integrals of energy, linear and angular momentum, and center-of-mass position. 
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